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Abstract 

Rare decays A& — > A7 and A5 — > Al + l~ (\— e, fi) are examined. We use QCD sum rules 
to calculate the hadronic matrix elements governing the decays. The A polarization in the 
decays is analyzed and it is shown that the polarization parameter in Aj, — > A7 does not 
depends on the values of hadronic form factors. And the energy spectrum of A in Aj, — > Al + l~ 
is given. 
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I. Introduction 

Processes associated with the flavor-changing neutral current (FCNC) b — ► s transition have regained much 
attention since the measurement of FCNC decays of the type b — > s 7 by CLEO [Q, |j . It is well known that these 
processes are forbidden at the tree level in the Standard Model(SM) and are strongly suppressed by the GIM 
mechanism, particularly for up type quarks. On the other hand. It is very sensitive to possible higher mass 
scales and interactions predicated by supersymmetic theories, two Higgs doublet models etc.. Such interactions 
shape the b — > s transition via operators and their Wilson coefficients appearing in the low energy AB = 1 
effective Hamiltonian. Hence the study of such process is an important way to test the CKM sector of the SM 
and possibly opens a window to physics beyond the SM. 

For the experimental side, there is some data on the exclusive decay B — > A'*7 and the inclusive decay 
B — > X s ^i with branching ratio Q 

Br(B+ -> J ftT+*(892) 7 ) = (5.7 ± 3.3) x 1CT 5 , 



Br(5° -> JT°*(892) 7 ) = (4.0 ± 1.9) x 1CT 5 , 

Br(5 -> X sl ) = (2.3 ± 0.7) x 10~ 4 . 

These data have prompted a number of study aimed at restricting the parameter space of various extensions of 
the SM Ji|. Similar analyses based on decays of B mesons have also been performed for the transition b — > sl + l~ , 
which has not been observed yet B and only upper limits on inclusive branching ratios have been given [Q : 



Br(6 -> se+e~) < 5.7 x 10" 5 , Br(6 -► sfi + fT) < 5.8 x 10~ 5 and Br(6 -► se ± /i =F ) < 2.2 x 10" 5 (at90%C*.L.). 

However, to analyze the helicity structure of the effective Hamiltonian mediating the transition b — > s, such 
analyses are not enough since the information on the handedness of the quark is lost in the hadronization 
process. To access the helicity of the quarks, analyzing the decay of baryons is the only way. One experimental 
drawback of baryon decay compared with B meson decay is that the production rate of A& baryons in b quark 
hadronization is (10.1 jl| '\)%, which is significantly less than that of B meson(i?r(& — * B + ) = Br(b — > B°) = 
(39. 7^2'®)%, Br(b — > B®) = (10.5^17)%) Hi > nence the relevant analyses have to wait for more data on heavy 
quark decays from future colliders. 

For exclusive heavy-to-light decays one has to calculate hadronic matrix elements of operators in the 
effective Hamiltonian between a heavy hadron and a light hadron, which is related to the nonperturbative 
aspect of QCD. There are a number of papers to calculate the hadronic matrix elements in exclusive rare decays 
B — > K(*H + l~ For exclusive heavy-to-light decays A& — * Al + l~ there are a lot of form factors to describe 

the hadronic matrix elements. However, for Ab we may use the heavy quark effective theory (HQET). It is well- 
known that HQET simplifies greatly the analysis of the decay of heavy hadrons in that heavy quark symmetries 
restrict the number of from factors, in particular, for the baryonic transition Aq — > light spin- 1/2 baryon, there 
are only two independent form factors irrelative to Dirac matrix of relevant operators p2j| ; for heavy hadron 
to heavy hadron transition, we only have one form factor, which is known as Isgur-Wise function |lO| . The 
computation of two form factors, F\ and F%, in HQET is the main work in analyzing exclusive decays of Ab to 
light baryon. Aj — > A7 has been in detail investigated by T. Mannel and S. Rocksiegel [Q], where the simple 
pole model was adopted to compute Fi and i*2- As we know, Ab — > Al + l~ has not been examined yet. The 
form factors F\ and F2 in A c — > A have been calculated in nonrelativistic and relativistic quark models 
In this paper we employ the widely applied approach of QCD sum rules, which is based on general features of 
QCD to calculate the F\ and F 2 . For our purpose, we can use some expressions of Ref. [jl3| due to the 
similarity. We analyze the A polarization. An interesting result is that the polarization parameter in Aj — * A7 
does not depends on the values of hadronic from factors, F\ and F 2 . 

This paper is organized as follows. In Sec. II we write down the SM effective Hamiltonian governing 
transitions b — > S7 and b — > sl + l~ , and give some information on the corresponding Wilson coefficients. In Sec. 
Ill, we compute the the Form factors F\ and F 2 , by using QCD sum rules. Sec. IV and Sec. V contribute 
to the analysis of decay Ab — > A7 and Ab — > Al + l~, respectively. In Sec. VI, we discuss our numerical results 
along with some relevant points. 
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II. Effective Hamiltonian 



In the SM, the Hamiltonian relevant for b — ► s transition consists of ten operators O^ i = 1, 10. Neglect- 
ing terms proportional to V u bV* s (the ratio \V u bV us * /V t bV ts *\ is of order 10~ 2 ), the effective Hamiltonian takes 
the form El 



where 
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Here a,/? are color indices, 6r.l = [(l±75)/2]6, and er A " / = (i/2)[7 A1 , 7"]); e and g s are the electronmagnetic 
and the strong coupling constants, respectively. The Wilson coefficients are given in Table I @|, where NDR 
denotes the naive dimensional regularization scheme and MS denotes the modified minimal subtraction scheme. 
The coefficients' dependence on regularization scheme must disappear in the decay amplitude if all corrections 
are taken into account. The operator responsible for decay b — > s 7 is O7, while operators responsible for decay 
b — > sl + l~ are Oi,O g ,O w - Because O7 does not include lepton fields, it has to be combined with 7ZZ vertex 
and lepton fields to contribute to b — > sl + l~. The result operator is O t = j^rribso M „(l + -f 5 )q" /q 2 bl'y tI i. 

We would like to write Og, Oiq and 7 in forms allowing for non-SM couplings (omitting the color indices): 
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(3) 



where hy,gv and h A ,g A are the vector and axial vector couplings respectively, in order to discuss possible 
effects of models beyond SM. In general the parameters hy and h A in Oiq may be different from those in 
Og and for the sake of simplicity we consider the case of same parameters as shown in Eq.(3). In the SM 
hv = l,h A = 1, fjv = 1,9a = —1 if we neglect the mass of strange quark in respect to b quark. 

As mentioned in Sec. I, the heavy quark symmetries restrict the number of form factors to two |12] |: 

< A(p,s)\sTb\A b (v,s) >= u A (p, s){F 1 {p.v)+ tfF 2 (p.v)}Tu Ab (v,s) (4) 
The following section contributes to the computation of F\ and F2 by QCD sum rules in HQET. 
III. Computation of F\ and F 2 

In Ref. |J, using QCD sum rules, we computed the form factors of to p transition. Here we follow the 
approach and use the some expressions obtained in the work. 
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To compute F\ and F 2 within the QCD sum rule approach we need to consider the three-point correlator 
U(P,P,z) = i 2 [ d 4 xd 4 ye lk - x -* p -y < 0\Tf(x)h v (0)Ts(0)j(y)\0 > , (5) 



of the flavor-changing current h v Ts and of A;, current j v and A current j, where P = m&« + k and z — P ■ v. 
The baryonic currents for Af, are 

f = e abc { q J a CTr q b 2 )h c v , 

where T has too kind of choice, 75 and 75 /b. We choose 75 for the sake of simplicity since the numerical 
differences resulting from the different choices of j v are not significant |1| . 



The three-quark tensor currents for proton and A are 17 



jT(xi,x 2) x 3 ) = a^ l5 d a (x 1 )u b (x 2 )Ca^u c (x 3 )e abc , 

Jl(xi,x 2 ,x 3 ) = a^ l5 d a (x 1 )u b (x 2 )Ca^s c (x 3 )e abc . 
After Fierz transformation, they can be written in S + P form 

J^(x 1 ,x 2 ,x 3 ) = 4[u a (x 3 )u b (x 2 )C 75 d c (x 1 )+ 75 u a (x 3 )u b (x 2 )Cd c (x 1 )]e abc , 

Jl(x u x 2 ,x 3 ) = 4[s a (x 3 )u b (x 2 )C l5 d c (x 1 )+ l5 s a (x 3 )u b (x 2 )Cd c (x 1 )}e abc . (6) 

The similarity allows us to use directly the analytical expressions obtained in the sum rule analysis of p final 
state to compute the form factors of A final state if we neglect the mass of s quark. 

After inserting a complete set of physical intermediate states, as the phenomenological consequence of (5), 
we have 

II(P' ,P,z) = f Ab fA , 2 t ; P+r[Fi (z) + F 2 (z) tf\ ^ + m \ +res, (7) 



(oj — 2 A) zw "P 2 



in 



where P + = (1+ ^)/2, A = m\ b — nib, uj — 2k ■ v and /a 6 , /a are the so-called "decay constants" which are 
defined by 

< Q\~f\A b >= f Ab u , < p(P)\j\0 >= f A u(P). (8) 

They can be found in Ref. |l5|, Q and Refs. |l7], respectively. To obtain (7) we have taken into account (4) 
and the heavy quark limit. 

By introducing the assumption of quark-hadron duality, the contribution from higher resonant and contin- 
uum states can be treated as 

res= / P^(W) 

J D > {v-u){s~P 2 ) y> 

The region D is characterized by one or two continuum thresholds s c . From the theoretical point of 
view, II(P , P, z) is combination of the perturbative contribution and the condensate contribution. In order to 
incorporate the above assumption, we should express the perturbative term in the form of dispersion relation 

nj, ert ( W ,P 2 ,z) = J dvds {v A ^ Z J p2y [i = 1,2), (10) 

where i = 1,2 denote the different terms associated with Fi and F 2 , respectively. In the standard way, we 
employ a double Borcl Transformation co — * M, P 2 — > T in order to suppress the higher excited state and 
continuum state contributions. Because the SU 3 violation effects arc small in QCD sum rule analyses of three 



point functions for mesons |20|, we expect the effects are probably even smaller for baryons. Therefore, we 



neglect the s quark mass in calculations and quote from Ref ]13| the analytical expressions of p l [y, s, z), i = 1, 2 
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and the condensate contributions. Thus, the resulting Borel transformed sum rules for Fi and F 2 can be written 

as 



and 



3^(« S GG)J T/4 (1 - M) e -4/3(l-4/3/T)/M 2 -8^/(TM) d/3i 
^(a s GG) - M) A e -4/3(l-4/3/T)/M 2 -8^/(TA/) d/3j 



Plert = ^^[-2z 3 a 3 -(-s + z(u + 2z)) 3 + 3z 2 (-s + z(u + 2z))a% 

p 2 pert = -^^[s-2z 2 + z{-v + CT)] 2 [vs + %z 3 -A t z 2 (-2v + a)-2z(-v 2 + hs + va% 



(11) 



(12) 



with a = \J — 4s + (v + 2z) 2 . It should be mentioned here that we only retain the condensates with di- 
mension lower than 7. In the numerical analysis, the "decay constants" and some other constants we used 
are|||||||: 

m Ab = 5.64GeV,TOA = 1.116GeV, f Ab = V0.0003GeV 3 

f A = 0.0208GeV 3 , A = 0.79GeV, (qq) ~ -(0.23 GeV) 3 , (13) 
(a s GG) ~ 0.04 GeV 4 . 

Again, owing to the small difference between p and A both in their mass and "decay constant" , the results 
here are similar to those of Ref. With the threshold v c = 3.5GeV, we can have a reasonably good window 
for F t , where 1.5GeV < — = M < 1.9GeV. The results are given in Fig.l and Fig. 2 respectively, where 
the different curves correspond to different choices of the Borel parameters. The comments made in Ref. [ |D| 
concerning about the sum rules, such as the condensates dominance property(for about 55%) , still hold here. 
One can see from figures 1 and 2 that R = F2/F1 = —0.42 at q 2 = q^ax which is almost the same as that of p 
final state. 

V. Decay A b -> A7 

As analyzed in Ref. 0], the long distance contribution to A b — > A7 is negligibly small, so we shall only 
consider the short distance contribution to the decay. The matrix element of the operator O7 between the initial 
and final state is 

(A(p, S ), 7 (k,e)\Or\A b (v,s')) = -l^mtiA&s^sa^igv- g Al5 )b\K b {v, s / ))( 7 (fc,e)|F^|0). (14) 

oZir 

In order to study the helicity of the final quark, we are concentrated on the decay rate of unpolarized Ab baryons 
into light baryons with a definite spin directions. We obtain the decay rate from Eqs.(4) and (14) as follows: 

r = ^(y t:Vtb f{^fmlml h {l - x 2 ){ g A±&[(l - 2x 2 + X *)\ Fl \ 2 + 2{x - 2x 3 + x 5 ) Fl F 2 

+ (x 2 - 2.x 4 + x 6 )\F 2 \ 2 } + gv9A(v.s)[(2x - 2x 3 )\F 1 \ 2 + 2(2x 2 - 2x 4 )F 1 F 2 + (2x 3 - 2x 5 )|F 2 | 2 ]} (15) 

where x = m A /m Ab . Our result is exactly the same as that of Ref. with real F± and F 2 , as it should be. 
In order to compare the result with experiment we rewrite the rate in terms of the polarization variables as 
defined in Ref. Q 

r = r [i + ap-s A ] (16) 

where p is the unit momentum vector of the A and S A is its spin vector. From expression (15), it is straight- 
forward to find that 

2x 2g v g A 



{l + x 2 ){g 2 v +g 2 A Y 

L/ 7 (_r K 



(KVtbfiT^) *mlmUl - x 2 ) 3 (F 1 + F 2 xf{g 2 v + g\ ). (18) 
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Eq. (17) is different from that in Ref. 0] and means that the polarization parameter a is independent of the 
form factors F-y and F 2 . The conclusion that a does not depend on the hadronic structure has also been obtained 
in Ref. pl| . However, in Ref. pl| F 2 — has been assumed so that the conclusion is trivially obtained and 
the factor 2x/(l+a; 2 ) in Eq. (17) is missed there. The QCD sum rules analyses give that F± — 0.50 ± 0.03 
and F 2 = -0.10 ± 0.03 at the point p° = (m Ab 2 + ra A 2 )/(2ra A J = 2.93GeV. Although this point is to the 
disadvantage of HQET application, which justifies Eq.(4), for the recoil at this point is the largest, we still 
assume that the heavy quark symmetries are applicable to some extent at this point. Taking x = 0.20, we have 

a = Q ^.J^_ (19) 

and 



r = (1.06 ± 0.16) x 10- 17 (^ + <?A)GeV, (20) 

where the uncertainty is rooted in the uncertainties of F\ and F 2 . In the SM, we take gy = 1, g A = —1, the 
corresponding branching ratio is 

BR(A b — > A7) = (3.7 ± 0.5) x 10~ 5 (21) 
which is within the range obtained in Ref. 
V. Decay A b -> Al + l~ 

Thc decay of lcpton final state is a little more involved than that of 7 final state in the integration over 
phase space. The relevant matrix element of the process are: 

(A(p,s),l (pi, si) , l + (pa, s 2 )\0 7 \A b (v, s')> = —^m b {A(p,s)\sa^ v {gv - gAl5)^b\A b (v, s')} 

I67H _ 1 

(r( P i, Sl )J + (P2,s2)|W|o) 

e 2 

(A(p,s),l ( Pl , Sl ),l + (p 2 ,s 2 )\0 9 \A b (v,s')} = — -j(A(p,s)|s7^(/i v - h A j 5 )b\A b (v,s')) 

OZ7T 

(r( Pl , Sl )j + ( P2 , S2 )\£ri\o) 

e 2 

(A(p,s),l (p 1 ,si),l + (p2,s 2 )\Oi \A b (v,s')) = ^i(A(p, s^sj^hy - h A j 5 )b\A b (v, s')) 

(l-( Pl , Sl )J + ( P2 , S2 )\£ 1 ^ 5 £\0) (22) 

To obtain from the above expressions the differential width with respect to the energy E of A is a matter of 
algebra, we obtain 

^ - A(y) + s ■ vB{y) (23) 



where y = E/m\ b . The expressions of A{y) and B(y) are given in the appendix for they are somewhat 

y 9 



lengthy. In the expressions, we have replaced Cg with Cg , which takes the main long distance contribution 
into consideration by following formalism p4B2, E3| 



2 



a" = Cs + (3Cl + cm...) (~) 



where 



if s < 4x 2 and 



, / x r 4, 2 8 16 x 2 4 Ux 2 . 4x 2 . Ax 2 . 1/2l 

h(x, s) = -[-In. 2 - - - - - + - yf — - 1(2 + -)arctan(- - 1)-V*] 



r 4, 2 8 16 x 2 2 f 4^\ 4a; 2 . n . 1 + Jl - Ax 2 / s . 

Hx, s) = -{-Inx 2 - — - — — + 2 + In V ' - wr } 

9 27 9 s 9 V s s 1-^1- 4x 2 /s 
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if s > Ax 2 , with x = m c /m& and s = q 2 /m 2 (q 2 = M? + ,_ ). As in Sec. V, we write the total width in the form of 
Eq.(16) 



r [l + ap-S A ]. (25) 



Our calculation yields 



T = (0.0045#1 + 0.0045 ffv + 0.0075g A h A + 1.64/4 - 0.0069g v h v + 1.64/4) ' ( L00 ± 0.24) x 10 _1 'GeV 
0.0037g v h A + 0.0037.g A9v - 0.037g A h v - \.73h A h v 

a = ^ ^ ^ o- x (1.00 ± 0.08), (26 

0.0045cr| + 0.0045g v + 0.0075g A h A + 1.6Ah 2 A - 0.0069g v h v + lMh 2 v v h y ' 

where the uncertainties are mainly rooted in the uncertainties of F\ and i*2. In the SM, we take hy = 1, h A = 
1,2V =l,9A= -1, then 2r' gi ves th e total decay width (6.57±1.58) x 10~ 17 GeV and a = -0.54±0.04 . The 

differential widths, A(y) and *^^- B{y) (the spin-dependent term of the energy spectrum in the rest frame 
of At), are given in Fig. 3 and Fig. 4 respectively. 

VI. Discussions 

In this paper we have analyzed some features of the rare decays At — > A7 and Af, — > Al + l~, using an 
approach based on three-point function QCD sum rules to compute the relevant form factors. 

We have considered only the short distance contribution for Af, — > A7. There are some estimations on 
the long distance contribution in Ref. 0], the results turn out to be negligibly small, the decay A& — > A7 is 
dominated by the short distance piece. For Af, — > Al + l~ we have taken into account the main long distance 
contributions included in the coefficient Cg . 

As we have mentioned in Sec. IV, in applying QCD sum rules analysis for A& — > A7, we need to assume 
the the heavy quark symmetry is applicable at its most disadvantage phase space corner. For the process 
Ab — > Art~, the heavy quark symmetry works well in most of phase space. In particular, one can see from the 
figure 1 that the form factor F\ becomes larger when z approaches the zero recoil point so that its contributions 
to the decay are dominant. 

Our results show that the total decay width T(Ab — > Al + l~) is larger than that for Ah — > 7, which is due to 
the dominance of the long distance (resonance) contributions to the decay At, — ► A/ + ^~. From the expressions of 
Eq.(19) and Eq.(26), in addition to that the total width of lepton final state is greater than that of 7 final state, 
its polarization effect is more remarkable than that of 7 final state either. So the measurements of polarization 
parameter in the decay A^ — > Al + l~ is of a complement to those in A^ — > A7 in order to discover new physics . 
From Eq.(26), we can also find that the contributions to decay Ah — > Al + l~ mainly come from O9 and 0\q. 

We emphasize that the polarization parameter in A;, — > A7 is independent of the hadronic structure of A 
in the heavy quark limit so that it is a good quantity to probe new physics beyond SM. In addition to non-SM 
couplings in 07, we have also considered the non-SM couplings in Og and O10 and examined their effects in rare 
decays of A&. It is expected that the measurements of A& — > A7 and Ah — > Al + l~ will give us some information 
on physics beyond the SM. 



[1] M. R. Ahmady, D. Liu and A. H. Fariborz, Z. Phys.C71, 301 (1996); T. M. Aliev, D. A. Demir, E. Iltan and N. K. 

Pak, Phys. Rev D54, 857 (1996). 
[2] T. Mannel and S. Recksiegel, J. Phys. G24, 979 (1998). 
[3] Particle Data Group, Eur. Phys. J. C3 , 1 (1998). 
[4] CLEO Collaboration, |hep-ex/9710003 . 

[5] D. Liu, Phys. Lett. B346, 355 (1995); Phys. Rev. D52, 5056 (1998). 

[6] P. Colangelo, F. De Fazio, P. Santorelli, and E. Scrimieri, Phys. Rev. D53, 3672 (1996), and references therein. 
[7] Hai-Yang Cheng and B. Tseng, Phys. Rev. D53, 1457 (1996). 

[8] M. A. Ivanov et al., Phys. Rev. D56, 348 (1997); R. Akhoury and W. Loinaz, Phys. Rev. D53, 1416 (1996). 
[9] For a review of the implications of the measurement of b — > S7 see J. H. Hewett, in Spin Structure in High Energy 
Processes, Proceedings of the 21st SLAC Summer Institute on Particle Physics, Stanford, California, edited by L. 



7 



1994),p. 463. Phys.Rev.Lett.71, 496 (1993); J. H. 



De Porc el and C. Dunwo odie (SLAC Report No. 444, Stanford, 
Hewett, |hep-ph/980337C| . 
N. Isgur and M. Wise, Nucl. Phys. B348, 276 (1991). 

M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl.Phys. B147, 385 (1979); B147, 448 (1979). 

T. Mannel, W. Roberts, Z. Ryzak, Nucl. Phys. B 355, 38 (1991). 

C. S. Huang, C. F. Qiao and H. G. Yan, Phys. Lett. B437, 40 (1998). 

B. Grinstein, M. J. Savage, and M. G. Wise, Nucl. Phys. B319, 271 (1989). 

Y. B. Dai, D. S. Huang, D. Liu and CD. Lii, Phys. Lett. B371, 99 (1996). 

A. Buras and M. Miinz, Phys. Rev. D 52, 186 (1995); A. Buras, |hep-ph/9806471 . 

E.V. Shuryak, Nucl. Phys. B 198, 83 (1982). 

W-Y. P. Hwang and Kwei-Chou Yang, Phys. Rev. D49, 460 (1994). 
V.A. Nesterenko and A. V. Radyushkin, Phys. Lett. B 128, 439 (1983). 
T. Huang and C. W. Luo, Phys. Rev. D 50, 5775 (1994L 



[23] 



hep-ph/9808431. 



Chun-Khiang Chua, Xiao-Gang He, and Wei-shu Hou, 
C. S. Lim, T. Morozumi, and A. I. Sanda, Phys. Lett. B 218, 343 (1989); N. G. Deshpande, and K. Panose. Phys. 
Rev. D 39, 1461 (1989. 

P. J. O 'Donnell and H. K. K. Tung, Phys. Rev. D 43, R2067 (1991); N. Paver and Riazuddin, ibid. 45, 978 (1992); 
A. Ali, T. Mannel and T. Morozumi, Phys. Lett. B 273, 505 (1991). 



8 



Appendix: Expression of A(y) and B(y) 



A(y) = h 2 m\ b {C 7 2 m b 2 r{F 1 2 [g A 2 {A - 8x 2 + Ax 4 ) + g v 2 (A - 8x 2 + Ax 4 )] + F 1 F 2 [g A 2 {8x - 16a; 3 + 8x 5 ) + 

g v 2 (8x - 16a; 3 + 8x 5 )} + F 2 2 [g A 2 (Ax 2 - 8x 4 + Ax 6 ) + g v 2 {Ax 2 - 8a; 4 + 4a; 6 )]) + C 1 2 m h 2 {F 1 2 [ 
g v {-A - Ylx - Ax 2 + Ay) + g A {-A + 12x - Ax 2 + Ay)} + FiF 2 [.g y 2 (16x - 8x 3 - 2Ay - 16xy) 
+g A 2 {16x - 8x 3 + 2Ay - 16xy)} + F 2 2 \g v 2 {-\2x + 12a; 2 - Ax 4 + I2y - 8x 2 y - 16y 2 ) + g A {l2x 
+12x 2 - Ax 4 + 12y - 8x 2 y - I6y 2 )}} + C 7 (Cg ff + Cg // *)m,,m Ab {Fi 2 [ 3A /i yl (-6x + 6x 2 - 6y + 6xy) 
+gvh v (~6x - 6a; 2 + 6y + 6xy)} + F 2 2 [g v h v (-6x + 6a; 2 + 6y + 6xy - 12y 2 ) + g A h A (-6x - 6x 2 
-6y + 6xy + 12y 2 )} + F^^gyhviUx - 12y - 12xy + 12y 2 ) - g A h A (12x + 12y - 12xy - 12y 2 )}} 
+ (Cio 2 + |Cg // | 2 )m 2 ib {Fi 2 [/i A 2 (3x - 2a; 2 + 3x 3 + 3y - 6xy + 3x 2 y - Ay 2 ) + h v 2 (-3x - 2x 2 - 3x 3 
+3y + 6xy + 3x 2 y - Ay 2 )} + FiF 2 [/ia 2 (6x + 2x 3 + 6y - 12xy + 6x 2 y - I2y 2 + 4a;y 2 ) + h v 2 (6x 
+2x 3 -6y- 12xy - 6x 2 y + 12y 2 + Axy 2 )} + F 2 2 [h A {3x + 2x 2 + 3x 3 + 3y - 6xy - x 2 y - 8y 2 
+Ay 3 ) + h v 2 (-3x + 2x 2 - 3x 3 + 3y + 6xy - x 2 y - 8y 2 + Ay 3 )}}} 

B(y) = -k 2 rxm\ b {-C 7 {Cl ff + Cl ff *)m b m Ab [6F 1 F 2 (g v h A - g A h v )x + F 2 (3g v h A - 3g A h v + g v h A x 
+g A h v x) + F 2 (-3g v h A + 3g A h v + g v h A x + g A h v x) + 2F 2 (F 1 g v h A + 3F 2 g v h A + Fxg A h v 
-2>F 2 g A h v )y} + (C w 2 + \C e 9 Si ' \ 2 )m%h A hv[-AF x F 2 x - F 2 (3 + x 2 ) + F 2 (3 + x 2 ) + 4(F 2 - 2F 2 
+F 1 F 2 x)y + AF 2 y 2 } - AC 7 2 t g A g v ml[-8F 1 F 2 x + F 2 (3 - x 2 ) + F 2 (-3 + x 2 ) + 2(F 2 - 5F 2 
+AF 1 F 2 x)y + 8F 2 y 2 )}} (27) 



where r = 1/(1 + x 2 - 2y), k 2 - \V bt V st \ 2 (aG F / ir) 2 ^y^x> / (I28n 3 ) . 

TABLE I. Wilson coefficients CAy) for Al!L = 225 MeV, y = 5 GeV and m t = 174 GeV. NDR scheme. 

V / MS 

Ci -0.243 

C" 2 1.105 

C 3 1.083 x 10 -2 

d -2.514 x 10~ 2 

C 5 7.266 x 10~ 3 

C 6 -3.063 x 10~ 2 

C* 7 -0.312 

C 9 4.193 

Cio -4.578 



FIG. 1. Form factor Fi(z) with different values of Borel parameter T and M. 



FIG. 2. Form factor F2(z) with different values of Borel parameter T and M. 



FIG. 3. The spin-independent term of the spectrum: A(y) 



FIG. 4. The spin-dependent term of the spectrum: ^ — —B(y) 
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